Effects of alloying and strain on the magnetic properties of Fei 6 N 2 
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The electronic structure and magnetic properties of pure and doped Fei6N2 systems have been 
studied in the local-density (LDA) and quasiparticle self-consistent GW approximations. The GW 
magnetic moment of pure FeigN2 is somewhat larger compared to LDA but not anomalously large. 
The effects of doping on magnetic moment and exchange coupling were analyzed using the coherent 
potential approximation. The theoretical Curie temperature in pure Fei6N2 is significantly higher 
than the measured value, which is attributed to the quality of available samples and the inter- 
pretation of experimental results. We found that different Fe sites contribute very differently to 
the magnetocrystalline anisotropy energy (MAE), which offers a way to increase MAE by small 
additions of Co or Ti. MAE also increases under tetragonal strain. 



I. INTRODUCTION 

Ordered nitrogen martensite a"-FeigN2 was first syn- 
thesized in bulk form by quenching of the cubic nitrogen 
austenite 7-FeN with a subsequent annealing.- Quench- 
ing initially produces disordered cv'-FeN, which then or- 
ders during low-temperature annealing to produce a"- 
Fei6N2. The latter is a metastable phase with a dis- 
torted body-centered tetragonal structure, which decom- 
poses into cv-Fe and Fe4N near 500 K. 

Interest in cv"-Fei6N2 was revived much later when 
it was synthesized in thin film form and a very large 
value (^3/i b) for the average Fe magnetic moment was 
reported. 2 This result was not independently confirmed 
until twenty years laterP Owing to the rapid develop- 
ment of the magnetic recording technologies, this con- 
firmation inspired numerous studies of thin-film samples. 
However, the existence of the "giant" Fe moment remains 
controversial as many researchers did not reproduce these 
findings, while others confirmed thempHSl The lack of 
consistent and reproducible experimental results may be 
attributed to the difficulties associated with the prepa- 
ration of single-crystal Fei@N2 and stabilization of ni- 
trogen, as well as with the accurate measurement of the 
magnetization in multi-phase Fe nitride samples. This is- 
sue has recently attracted additional interest due to the 
search for new permanent magnetic materials without 
rare-earth elements!^ A new way to prepare single-phase 
FeigN 2 powder was recently reported, along with evi- 
dence of high maximum energy product (BH max )W 

Most theoretical studies of the magnetization of a"- 
Fei6N2 were performed using the local density approx- 
imation, generalized gradient approximation (GGA) or 
LDA+U, though recently Sims et al^ applied a hybrid 
functional and the GW approximation to this material. 
In LDA or GGA the magnetic moment of FC16N2 is 
only slightly enhanced compared to elemental Fe. Lai 
et ati^l included electronic correlations within LDA+Z7 
and found an enhanced magnetization Af=2.85/is/Fe. 



Wang et aZp^KQl identified a localized Fe state coexist- 
ing with the itinerant states in X-ray magnetic circular 
dichroism (XMCD) measurements. They introduced a 
specific charge transfer between different Fe sites and ob- 
tained a large M in LDA+U. However, the choice of the 
correlated orbitals and the associated value of the Hub- 
bard U parameter is not well-defined for metallic systems. 
For example, the on-site interaction parameters obtained 
by Sims et al. 9 using the constrained RPA differ substan- 
tially from those proposed by Wang et al.. The quasipar- 
ticle self-consistent GW approximation (QSGW) 14 * 15 ! is 
more reliable and provides a more satisfactory way to de- 
termine the ground state density and magnetic moment. 
In the present paper we apply this method to Fei6N2. 

Studies of exchange interaction, Curie temperature 
(T c ), and MAE of Fei 6 N 2 met with additional difficul- 
ties. In particular, measurements of Tc are hampered by 
the decomposition of the metastable FeigN 2 into Fe4N 
and Fe, which was reported to occur above 200° CP in 
the 230-300°C rangepor at 400°CPSugita et al. extrap- 
olated their data to estimate Tq at 540°CP Thermal sta- 
bility of Fe_i6N 2 was reported to increase with addition of 
Co and TpEl ( up to 700°C in the Ti case). However, no 
experimental information is presently available about the 
Tc of Co or Ti-doped Fei6N 2 , or of any other Fei6N 2 sam- 
ples stabilized at high temperatures. To the best of our 
knowledge, there have been no theoretical studies of the 
exchange interaction and Curie temperature in FeigN2. 
Systematic studies of the effects of doping on M and Tq 
in Fei6N2 also appear to be lacking. 

As for the MAE, only a few experimental values were 
reported, and they are varied and inconclusive. For ex- 
ample, Sugita et alW obtained an in-plane MAE, while 
Takahashi 18 found a large uniaxial MAE. The only avail- 
able theoretical calculations of MAE used an empirical 
tight binding (TB) approximation! 1 ^ 

In this paper we study the magnetization, Curie tem- 
perature, and magnetocrystalline anisotropy energy of 
pure and doped Fei6N 2 using several well-tested elec- 



tronic structure techniques and suggest possible routes 
for improving its properties for permanent magnet appli- 
cations. 



II. COMPUTATIONAL METHODS 

Most LDA, GGA and QSGW calculations were per- 
formed using a full-potential generalization^ of the stan- 
dard linear muffin-tin orbital (LMTO) basis seti^H This 
scheme employs generalized Hankel functions as the enve- 
lope functions. Calculations of MAE we also performed 
using the recently-developed mixed-basis full-potential 
method,^ which employs a combination of augmented 
plane waves and generalized muffin-tin orbitals to repre- 
sent the wave functions. The results of a traditional non- 
self-consistent application of the GW approximation de- 
pend on the non-interacting Hamiltonian generating the 
self-energy. This issue can be particularly problematic 
for metals. In contrast, QSGW method does not suffer 
from this limitation: it is more reliable than the standard 
GW. This method gives quasiparticle energies, spin mo- 
ments, dielectric functions, and a host of other properties 
in good agreement with experiments for a wide range of 
materials, including correlated ones such as NiO. The de- 
tails of QSGW implementatiorfHESl and applications can 
be found elsewhere. 

The pair exchange parameters were obtained using two 
linear response approaches: 

(1) Static linear-response approaches implemented 
within the atomic sphere approximation (ASA) to the 
Green's function (GF) LMTO methodP In addition to 
making a spherical approximation for the potential, this 
method makes the long- wave approximation (LWA), so 
that the pair exchange parameter is proportional to the 
corresponding spin susceptibility Xij ■ The exchange 
parameters Aij obtained in this method are related to 
the parameters of the classical Heisenberg model 



H — / y Jij Oj ' bj , 



(1) 



by the following renormalization for ferromagnetic(FM) 
and antiferromagnetic(AFM) cases: 



= -4Aij/tOimj = 



(2) 



AAij/miirij FM 
-AAij/rmmj AFM 



where m^ is the magnetic moment on site i. With this 
renormalization all results obtained for the Heisenberg 
model Eq.(fT|) can be used directly. 

(2) Dynamical linear response approach with the bare 
susceptibility x(q, w) calculated in the full product basis 
set representation using the LDA or QSGW electronic 
structurepH The results are then projected onto the 
functions representing local spin densities on each mag- 
netic site, which gives a matrix Xy(q, w) in basis site 



indices.-- 1 This projection corresponds to the rigid spin 
approximation. The inversion of this matrix with a sub- 
sequent Fourier transform provides the real-space rep- 
resentation of the inverse susceptibility representing the 
effective pair exchange parameters: 

Jii = lim n 7T^ / <*q Wq, ^P 1 e' qR - • (3) 

w — >0il BZ J 

Tc is calculated both in the mean-fielcP^ (MFA) and 
the RPaISS approximations. The actual Tq may usually 
be expected to lie between the results of these two ap- 
proximations. 

To address the effects of doping, we used our im- 
plementation of the coherent potential approximation 
(CPA) within the TB-LMTO code, which follows the 
formulation of Turek et a/P" and Kudrnovsky et al.W^ 
A coherent interactor matrix fli is introduced for each 
basis site i treated within CPA. At self-consistency ga — 
(Tj — Qj)~ , where 7i is the coherent potential matrix for 
site i, and ga is the on-site block of the average auxiliary 
LMTO GF matrix g = (? - S)' 1 . This on-site block is 
extracted from the Brillouin zone integral of g(k). The 
conditionally averaged GF at site i occupied by compo- 
nent a is gfi = (P a — fi) _1 , and the CPA self-consistency 
condition can be written as ga = Yl a c i9iii here cf is 
the concentration of component a at site i. Using this 
equation, at the beginning of each iteration the stored 
matrices fii are used to obtain an initial approximation 
to CP,;. In turn, 3^ is used in the calculation of ga by a 
Brillouin zone integral. The next approximation for fli 
is obtained from Hi = 3^ — g^ . These output matrices 
are then linearly mixed with the input fli matrices at the 
end of the iteration. 

We found that the mixing coefficient of 0.4 for $7^ works 
well in most cases. For fastest overall convergence, we 
found that it is usually desirable to iterate CPA iterations 
until the f^ matrices are converged to a small tolerance, 
and only then perform the charge iteration. The conver- 
gence of $1 is done separately for each point on the com- 
plex contour to the same tolerance. With this procedure, 
fairly aggressive Broyden mixing can be used for LMTO 
charge moments. CPA convergence at each charge it- 
eration usually takes 10-50 iterations depending on the 
imaginary part of energy and the selected tolerance. At 
the beginning of the calculation, the f2j matrices are set 
to zero; afterwards they are stored and reused for subse- 
quent iterations. In order to avoid unphysical symmetry- 
breaking CPA solutions (which otherwise often appear), 
the coherent potentials and the k-integrated average aux- 
iliary GF are explicitly symmetrized using the full space 
group of the crystal. As a result, the use of CPA does not 
impose any restrictions on the symmetry of the crystal. 
Calculations reported here were performed without using 
charge screening corrections for the Madelung potentials 
and total energy. 

The effective exchange coupling in CPA is calculated 
as 



J (c) = cJ A (c) + (1 - c)J B (c) 



(4) 



where the component-specific Ji (c) are obtained using 
the conditionally averaged GF and the formalism of Ref. 

ED 

For MAE calculations the self-consistent solutions are 
found including spin-orbit coupling (SOC) terms of or- 
der 1/c 2 . The MAE is defined below as K = E 100 ~E 001 , 
where -Eooi an d -E100 are the total energies for the mag- 
netization oriented along the [OOf] and [100] directions, 
respectively. Positive (negative) K corresponds to uni- 
axial (planar) anisotropy. We used a 24 x 24 x 24 fc-point 
mesh for MAE calculations to ensure sufficient conver- 
gence; MAE changed by less than 2% when a denser 
32 x 32 x 32 mesh was employed. All calculations except 
QSGW were performed with both LDA^ and GGA^ 
exchange-correlation potentials for comparison. 



III. RESULTS AND DISCUSSION 

The crystal structure of FeigN 2 is body-center- 
tetragonal (bet) with space group I 4 / mmm (#139). It may 
be viewed as a distorted 2x2x2 bct-Fe superlattice with 
c/a— 1.1. Crystal structure of FeigN2 was first identi- 
fied by Jack.W Here we use lattice constants a=b=5.72A 
c=6.29A and atomic position parameters Z4 e =0.3125 and 
£8^=0.25 from Jack's work as the experimental structure, 
(see Fig. [I]). 

We also relaxed the structure by minimizing the total 
energy in LDA and obtained Z4 e =0.293 and :r8/i=0.242, 
nearly identical to that obtained by Sawada et apH. The 
primitive cell contains one N and eight Fe atoms divided 
into three groups indicated by Wyckoff sites: two 4e, 
four 8h and two Ad sites (correspondingly first, second 
and third neighbors to N). 



A. Magnetic moments and electronic structure 

Table IT] shows the atomic spin moment rrii at the 
three Fe sites and magnetization M (orbital magnetic mo- 
ment is small, hereafter we only include spin magne- 
tization in M). Within the LDA, M=2.38ju B /Fe was 
obtained, in good agreement with previously reported 
calculations 4 " 6 . The enhancement relative to elemental 
bcc-Fe has been attributed to the size effeclpH QSGW 
gives M=2.59fj, B /Fe, about 9% larger than LDA. While 
it is known that GW enhances spin moments relative 
to LDA, to ensure the genuineness of this enhance- 
ment of moment, we also carried out the QSGW cal- 
culation of bcc-Fe and found that QSGW enhance the 
LDA magnetic moment in elemental Fe by only ~2% 
(2.20— >2.24/is). Sims et al. 9 found a similar M in their 
GW calculation while they also obtained a larger mag- 
netization enhancement in bcc-Fe with M=2.65//s/Fe. 
Considering Fei6N 2 consists of about 87% Fe, the 9% en- 
hancement we find non-trivial. However, it is still well be- 
low M= 2.85^ B /Fe, obtained in LDA+C/ by Lai et aZPl 
The spin moment on the 4c? site reaches rrii= 3.11/is in 
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FIG. 1: (Color online) (a) Crystal structure of Fei6N2. The 
experimental atomic positions are shown. Relaxed structure 
have slightly different Z4 e and xsh- (b) {001} plane with 8/1 
and N atoms, (c) {110} plane with 4e, 8h and N atoms, (d) 
{110} plane with Ad and 8/1 atoms. 



TABLE I: Atomic spin magnetic moment mi and spin magne- 
tization M in Fei6N2 in different methods. Calculations are 
in the LDA unless GGA or QSGW is specified. 



Mpthod 


mi(fi B ) a 




M c 




4e 


8/1 


Ad 


N 


(hb/Fc) 


(emu/g) 


ASA 


2.07 


2.40 


3.03 


-0.06 


2.48 


2.47 


239 


ASA-GF 


2.10 


2.41 


2.99 


-0.10 


2.48 


2.47 


239 


FP 


2.08 


2.32 


2.84 


-0.05 


2.39 


2.38 


231 


FP(GGA) 


2.21 


2.40 


2.86 


-0.04 


2.47 


2.43 


236 


QSGW 


2.24 


2.55 


3.12 


-0.01 


2.62 


2.59 


251 



"Spin moment inside atomic or muffin-tin sphere. 

'Average of the atomic spin moments of all Fe sites without taking 
account of interstitial and N sites. 

c Average spin moments within the cell (with taking account of 
interstitial and N sites). 



QSGW, though we do not observe any obvious charge 
transfer from Ad to 4e and 8/1 sites in QSGW, relative to 
the LDA. Hence, we can not attribute the enhancement 
of M to the charge transfer as suggested by others^"^. 
Density of states (DOS) calculated within LDA and 
QSGW are shown in Fig. [2j The LDA result is simi- 
lar to previously reported results. A careful examination 
of the band structure reveals that QSGW significantly 
modifies the energy bands near Ep, relative to LDA. It 
has a slightly larger on-site exchange, widening the split 
between the majority and minority DOS and increasing 
M by about 9%. Both DOS figures show hybridization 
between N-2p and Fe-3d states at around —7 eV, indi- 
cating that QSGW does not strongly modify the relative 




FIG. 2: (Color online) Site and spin-projected densities of 
states within LDA (a) and QSGW (b). 



alignment of N-2p and Fe-3d levels. Comparing the par- 
tial DOS reveals that bands are slightly wider and hy- 
bridization is overestimated in LDA, as is typical since 
LDA tends to overestimate 3d bandwidths slightly. The 
DOS also show hybridization is stronger in the 4e and 
8h channels while weaker in the id channels, which are 
the furthest removed from N. Also, as typical with sec- 
ond row elements, QSGW pushes the N-2s bands down 
relative to LDA, from -16.2 eV to -18 eV. The N-2s also 
hybridizes with Fe-4e and Fe-8/i. However there is almost 
no hybridization with the furthest Fe-4<i at all because 
the N-2s orbitals is very localized. 



B. Exchange coupling and Curie temperature 

The Heisenberg model parameters Jij using LDA- ASA 
in the LWA, FP-LDA and FP-QSGW are plotted in 
Fig. [3] and tabulated in Table [TTJ The two LDA results are 
quite similar, confirming that the ASA and the LWA form 
a reasonable approximation. QSGW shows some differ- 
ences, particularly reducing those interactions which are 
AFM. 

The structure of J^ is much more complicated in 
Fei6N2 than in elemental bcc-Fe. The vectors connected 
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FIG. 3: (Color online) Real-space magnetic exchange param- 
eters J^ in Fei 6 N 2 within ASA-GF (a), FP (b), and QSGW 
(c) as functions of distance Rtj/a. (d) (Rij/a) 28 Jij in ASA- 
GF as a functions of distance Rij/a. The in-plane lattice 
constant a in Fei6N2 is as twice large as in bcc-Fe. 



the nearest 8/i-4e or 8h-4d sites are nearly along [111] 
direction, the magnetic interactions between them are 
generally large. In comparison, the largest interaction is 
also between the nearest sites connected by vectors along 
the [111] direction in bcc-Fe. Jij between Sh-Ae sites is 
very anisotropic due to the distortion of lattice around 
N atom. Similar anisotropy was also found for the in- 



TABLE II: Pairwise exchange parameters JV,(meV) and 
Tc(K) calculated with different methods. 





\R l3 \/a 


direction 


ASA 


FP 


GW 


4e-4e 


0.454 








-0.455 


19.37 


23.75 


16.75 




0.645 








0.645 


4.18 


14.20 


3.76 




0.713 


-0.5 


-0.5 


0.095 


3.24 


2.33 


1.88 




0.895 


-0.5 


-0.5 


-0.550 


-1.92 


-2.91 


-1.44 


8/1-8/1 


0.483 





-0.484 





2.00 


-2.55 


1.69 




0.516 





0.516 





1.49 


1.37 


4.84 




0.550 


0.016 


0.016 


0.550 


5.74 


4.60 


4.67 




0.684 


-0.484 


-0.484 





7.00 


11.07 


11.80 




0.707 


0.516 


-0.484 





0.20 


-0.68 


0.13 




0.730 


0.516 


0.516 





0.39 


-0.01 


-2.50 


4d-4d 


0.550 








0.550 


2.76 


1.54 


3.67 




0.707 


-0.5 


0.5 





1.15 


1.44 


0.78 




0.895 


-0.5 


-0.5 


0.550 


-2.37 


-4.49 


-1.90 


4e-8/i 


0.430 


-0.258 


-0.258 


-0.227 


20.55 


20.00 


14.08 




0.470 


0.242 


0.242 


0.323 


4.20 


5.88 


5.88 


8/i-4d 


0.448 


0.258 


-0.242 


0.275 


15.73 


22.30 


16.18 




0.827 


-0.242 


-0.742 


-0.275 


0.58 


0.60 


0.64 


4d-4e 


0.502 


-0.5 





0.048 


1.56 


1.85 


4.72 




0.708 





-0.5 


-0.502 


0.49 


-0.23 


-0.16 


T C (MFA 


) 






1552 


1621 


1840 


T C (RPA 


) 






1118 


1065 


1374 



plane couplings on the 8h lattice. Interestingly, a large 
coupling, (Jij— 23.75 meV in the FP-LDA calculation), 
occurs between two 4e atoms along [001]. This pair has 
been squeezed together by neighboring N atoms. Since 
exchange coupling is sensitive to the distance between 
those two sites, we also examined this exchange param- 
eter for the experimental atom coordinates, for which 
the bond length of the 4e-4e pair shrinks from 2.60 to 
2.36 A, and found that this Jij increases from 23.75 to 
36.55 meV, indicating significant exchange-striction ef- 
fect. The second nearest 4e-4e coupling(two Fe atoms 
with a N atom between them along (110) direction) is 
14.2 meV in FP and 4.18 meV in ASA. The relatively 
large disagreement may be a consequence of the shape 
approximation used in ASA, considering the presence of 
N atom and strong lattice distortion around this pair of 
atoms. 

The calculated magnetic interactions between different 
types of atoms have very different spatial dependence 
and correspondent asymptotic behavior. To demon- 
strate it explicitly on Fig. [3td) we show Jy scaled with 
(Rij/a) 2,s . With this renormalization J^ between atoms 
on 4e positions (smaller moments) are approximately 
constant in this range of distances(long-ranged interac- 
tion), while Jij between Fe atoms on 4d sites (with the 
largest moments) decay much faster( short-ranged inter- 
action), corresponding to more localized moment behav- 
ior. Such very different asymptotic behaviour suggests 
that these localized and delocalized interactions corre- 
spond to Fermi surface shapes with different dimension- 
alities. 

Tq is calculated from the exchange parameters and 
tabulated in Table [ill RPA values are about 30% smaller 



than the MFA ones. Typically experimental values fall 
between the MFA and RPA results, with the RPA being 
closer to the experiment in normal three dimensional sys- 
tems. In the present case, however, the reported extrap- 
olated experimental estimate Tc=810K-2- is smaller than 
both of MFA and RPA values, and smaller than the one in 
bcc-Fe^. Some portion of discrepancy may be attributed 
to the experinrental interpretation of the measured Tq. 
As noted in the introduction, Fei6N 2 decomposes with 
increasing T; moreover, there is a transition to the (non- 
magnetic) 7 phase at 1185 K. Since the measured M(T) 
is not a measurement of the single-phase material, it is 
still unknown what is the experimental value for the Tq 
in a single-phase Fei6N2. Unfortunately, disentangling 
the structural and magnetic degrees of freedom is very 
difficult, experimentally. Considering that M increases 
in Fei6N2, we expect Tq to increase also. 

Let us now discuss the influence on Tq when other 
atoms substitute for Fe or N. 

The CPA is an elegant, single-site ab initio approach 
to study substitutional alloys. We have implemented the 
CPA, including a MFA estimate for Tq, within the ASA. 
As we have seen by comparing exchange interactions in 
the ASA to those without this approximation, the ASA 
does not seem to be a serious approximation to the LDA 
in this material. 

Fig. [4] shows the M and and the normalized effec- 
tive exchange (or MFA estimation of Tq in units of pure 
Fei6N2) with doping by different elements. Both Co and 
Mn doping cause Tq to decrease. On the other hand, 
with N site being doped with B,C,P or Al elements, Tq 
and M change slightly and the Fermi surface character 
is barely affected. 

As shown in Fig. [4j Co or Mn-doped Fei@N2 decrease 
moment and exchange coupling. We neglected the pos- 
sible site preference effect in this calculation, and doped 
all three Fe sites with equal probability. Table |HI| shows 
the magnetic moment and Jo parameters of Fe and sub- 
stitutional components on all three different Fe sites. It 
indicates an opportunity to increase Tq by using a sepa- 
rate Co-doping on Fe-4d sites. 

Magnetic moments of the Fe component decrease with 
Mn doping and slightly increase with Co doping. With 
Co doping, the Fe moments on 8/1 and 4d sites do slightly 
increase, however this increment is not big enough to 
overcome the decrease resulting from Co substituting for 
Fe - the system behaves more like localized moments sys- 
tem. We also carried out the FP calculation of Fe7CoN, 
with one out of eight Fe atoms being replaced by Co 
atom and confirmed that the magnetization decrease, es- 
pecially when Co replace the Fe on 4e site. This is con- 
sistent with the CPA results. Mn substituent have larger 
magnetic moments than Co substituent. However, Mn 
doping decreases moments on Fe sites. Overall, the de- 
pendence of the total magnetic moments on substituent 
concentration are almost exactly same with Co and Mn 
doping. Another interesting observation is that magnetic 
moments of both substituents decrease with increasing of 



TABLE III: Component-resolved atomic spin moments mi ( the atomic spin moment of substitutional component are given in 
parentheses), magnetization M and exchanges Jo in Co and Mn-doped FeigN2 calculated within ASA-GF. 



Substituent 


X 


mi(fi B ) 


M(^s/atom) 


Jo(meV) 


4e 


8ft 


Ad 


4e 


8ft 


Ad 


Co 


0.00 
0.10 
0.20 


2.10 ( 1.44) 
2.09 ( 1.34) 
2.08 ( 1.27) 


2.41 ( 1.68) 
2.43 ( 1.64) 
2.45 ( 1.62) 


2.99 ( 2.11) 

3.01 ( 2.09) 

3.02 ( 2.08) 


2.47 
2.40 
2.32 


12.95 (11.37) 
11.97 ( 9.62) 
11.24 ( 8.47) 


15.70 (14.62) 
15.45 (13.56) 
15.19 (12.91) 


16.96 (19.42) 
16.70 (18.68) 
16.50 (18.07) 


Mn 


0.00 
0.10 
0.20 


2.10 ( 1.88) 
2.07 ( 1.75) 
2.06 ( 1.64) 


2.41 ( 2.25) 
2.35 ( 2.04) 
2.30 ( 1.85) 
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12.96 ( 7.05) 
12.03 ( 5.64) 
11.15 ( 4.41) 
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FIG. 4: (Color online) Spin magnetization M (a) 

and normalized exchange Jo/Jo(FeigN2) (with respect to 
pure FeieN2) (b) as functions of doping concentration in 
Fei6N2. The concentration x of doping element T is de- 
fined as (Fei_ a! T :E )i6N2 with Fe site doping (T=Co,Mn); and 
Fei 6 (Ni_ :E T a; )2 with N site doping(T=B,C,P and Al) . 



doping concentration. This can be explained by the par- 
tial density of states as shown in Fig. [5] The magnetic 
moment of Co slightly decreases with increasing of dop- 
ing concentration. As shown in Fig. [5] the unoccupied 
DOS peak right above the Fermi energy (E-p) in the mi- 
nority spin channel moves toward it. More electrons fill 
in the minority channel and decrease the magnetic mo- 
ment as doping increases. With Mn doping on the other 
hand, the peak in the majority channel right below E-p 
becomes less pronounced. It shifts toward Ep and de- 
creases the magnetic moment of Mn component. For the 
Fe component DOS, there is no peak structure near the 
Ep, and magnetic moment is much less sensitive to the 
substitutional concentration. 



FIG. 5: (Color online) Densities of state of substitu- 

tional component in random alloy (Fei-^Co^ie^ (a) and 
(Fei-zMn^ieNa (b). 



TABLE IV: Previous works on magnetic anisotropy in Fei6N2. 



Method 




A'(i^F) 


Easy axis 


Ref. 


Exp. 


Sugita et al. 


4.8 


[100] 


La] 




Takahashi et al. 


200° 


[001] 


18 




Takahashi et al. 


97 


[001] 


19 




Kita et al. 


44 


[001] 


36 




Ji et al. 


100 6 


[001] 


37 


TB c 


Uchida et al. 


140 


[001] 


|19| 



"Value of {K±+K 2 ). 

'Measured in partial-ordering FeigN2, author claimed MAE 
should be much higher for the single-phase sample. 
c Tight binding approximation 



C. Magnetic anisotropy 

Values of MAE from previous work are summarized 
in Tables |IV| Results of present work are shown in Ta- 
ble [V] All calculations are carried out within LDA unless 
GGA is specified. For the pure FeigN2, both experi- 



TABLE V: The MAE K, on-site orbital magnetic moment I 
and the AMAE A with different spin quantization axis direc- 
tion in pure, Co-doped and Ti-doped Fei6N2. Spin quantiza- 
tion axis direction e are along [001], [100] and [110] directions 
respectively. With the spin along [100] and [110], A and K 
values(with respect to [001] direction) directions are given. 
To estimate A, ^=50,70 meV had been used for Fe and Co 
atoms respectively. 





e 


K 


/(10-^s) 


A 




ficV 
Pc 


10" org 
era' 1 


4e 


8h 


Ad 


MeV 

Fc 


Exp." 


001 






54 


45 


71 






100 


116 


144 


35 


49 


64 






110 


116 


144 


36 


58 39 


64 




Exp." 


001 






52 


44 


68 




GGA 


100 


105 


131 


36 


48 


61 


110 




110 


105 


131 


36 


57 39 


61 


110 


Theo." 


001 






62 


46 


67 






100 


84 


103 


39 


50 


63 


137 




110 


84 


103 


39 


58 41 


63 


137 


Theo." 


001 






56 


43 


62 




GGA 


100 


52 


65 


38 


47 


58 






110 


52 


65 


38 


55 40 


58 




Fe 7 CoN 


001 






91 d 72 


47 


76 




(4e) c 


100 


165 


206 


49 d 30 


49 


61 64 


337 




110 


165 


206 


49 d 30 


61 36 


63 


336 




001 






63 


69 d 41 41 44 


70 




(8h) 


100 


42 


52 


36 


77 d 46 50 47 


67 


123 




110 


16 


20 


40 


90 d 38 38 56 


67 


81 




001 






63 


51 


120 d 81 




(Ad) 


100 


138 


171 


33 38 


52 


106 d 70 


271 




110 


138 


171 


36 


62 42 


106 d 70 


271 


Fe 7 TiN 


001 






ll d 63 


40 


69 




(4e) 


100 


127 


158 


10 d 27 


46 


65 65 


62 




110 


127 


158 


10 d 27 


55 38 


65 


62 




001 






55 


14 d 48 48 41 


69 




(8h) 


100 


57 


71 


38 


13 d 48 50 50 


65 


122 




110 


43 


53 


35 


13 d 41 41 62 


68 


83 




001 






60 


40 


14 d 69 




(Ad) 


100 


102 


127 


38 38 


45 


13 d 67 


83 




110 


103 


128 


38 


52 38 


13 d 67 


83 



"Exp. the experimental crystal structure was used. 

^Theo. the theoretically optimized crystal structure was used. 

c Doping site of the substitutional atom. 

d Orbital magnetic moments of the substitutional atoms. 



mental and optimized structure are investigated. Note 
that doping and SOC lower the symmetry, and the de- 
generacy of I varied. Within LDA, a uniaxial magnetic 
anisotropy K=14Ax 10 5 erg/cm 3 was obtained with ex- 
perimental atomic coordinates. Structural optimization 
gives a smaller MAE with K= 103xl0 5 erg/cm 3 . GGA 
gives smaller MAE than LDA. It is usually non-trivial 
to analyze the origin or site dependence of magnetic 
anisotropy. Fortunately in this system the value of spin 
orbital coupling on Fe atoms (^50 meV) is much smaller 
than exchange splitting (~2 eV) and one can safely use 
perturbation theory analysis^ where MAE is associated 
with anisotropy of orbital magnetic moment. In this 
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FIG. 6: (Color online) AMAE on 4e, 8ft, Ad Fe sites (A;) and 
their average value (A) and the LDA total energy (£"lda) as 
functions of spin quantization axis rotation. 



approximation the atomic magnetic anisotropy energy 
(AMAE) Aj is defined as the difference of the SOC en- 
ergies along different magnetic field directions, that is in 
turn defined by the corresponding anisotropy of orbital 
magnetic moments: 

k ~ J2 Me = 90°) = J2 &Mii 01 - ^ 100 )/4 (5) 

where £j is a SOC parameter, while m, and li are 
atomic spin and orbital magnetic moments correspond- 
ingly. Sum of Aj can be compared with total MAE K 
obtained using the total energies. This perturbation esti- 
mate takes into account of the fact that the spin moment 
has very weak anisotropy and the main change in £L • S 
product comes from the change of orbital magnetic mo- 
ment. This is the case for FeigN2 and we show some 
orbital magnetic moments I in Table |VJ For the pure 
Fei6N2, when spin quantization axis rotates from [100] 
to [001], I decreases on 8h sites, but increases on Ad and 
Ae. While the I on different sites behave differently, the 
total I increases during this rotation, which agrees with 
the predicted uniaxial character of MAE. When the spin 
quantization axis is along [110], SOC lowers the symme- 
try, and splits the four equivalent 8h sites into two pairs 
with I increasing on one pair and decreasing on the other. 
As shown in Fig. [61 there is a strong correlation be- 



tween K and A (with respect to magnetic field along 
[001] direction, and atomic value £i=50meV is used for all 
three different Fe sites for simplicity), where i indicates 
all atomic sites. Obviously, the atomic 8h sites make 
negative contributions to the desired uniaxial MAE, and 
while Ac and Ad sites make positive contributions. Thus, 
one may hope that doping on 8h site, thus eliminating 
negative (in-plane) contribution to MAE, may improve 
the uniaxial MAE. 

Since Co and Ti do ping had been reported to stabilize 
the F816N2 phase^E^, it seems logical to study prediction 
above using these dopants. We replaced one out of eight 
Fe atoms in the primitive cell with Co or Ti atom and 
relax the atomic positions within LDA and then study 
the anisotropy. If we replace one of four 8h atoms with 
Co atom, we found that Co atom has a larger I than any 
other Fe sites, however, it does not eliminate the negative 
contribution from 8h sites. Instead, it makes K smaller. 
Also I and then K along [100] and [110] directions become 
more anisotropic. Surprisingly, however, with a Co atom 
on 4e or 8h sites, the I difference between out-of-plane 
and in-plane cases become even larger on Ac and Ad sites 
and smaller on 8h sites. In other words, it makes the 
positive contribution from 4e and Ad sites stronger and 
the negative contribution from 8h sites smaller. As a re- 
sult, calculated MAE is doubled ( if=206xl0 5 erg/cm 3 ) 
within LDA with doped Co being on 4e site. A similar 
effect had been found with Ti doping. K increases when 
Ti is substituted on the 4e and Ad sites and decreases 
when substituted on 8h. Unlike the Co doping, magnetic 
orbital moments of Ti atom are small and barely change 
with spin rotation. Generally, for the same structure, K 
is always strongly correlated with A. The larger A is, 
the larger K is along that specific spin quantization di- 
rection. However, this correlation may not longer hold 
true with different structures. For example, in Fe7TiN, 
A is the largest with Ti doped on 8h site, however K 
is much smaller than those with Ti doped on 4e and Ad 
sites. 

Tctragonality is another factor which may affect the 
anisotropy in a significant way. Let us compare Fe 16 N 2 
with bct-Fe, where even for c/a=l.l (the c/a ratio for 
Fei 6 N 2 ) MAE is still rather tuxyP, In Fig. [7] the cal- 
culated MAE in Fei6N2 is shown as a function of c/a. 
This dependence is much stronger than in bct-Fe and 
we assume that MAE mostly originates from distortion 
of Fe sublattice around the N atom and the Fe-N hy- 
bridization. Experimentally, the large tetragonality can 
be obtained in films, where it can be tuned by the nitro- 
gen concentration 37 . However, according to our results 
above, doping bulk Fei6N2 in a way that increases c/a 
may lead to MAE increase. The MAE and AMAE are 
well correlated as shown in Fig.[7| Within this c/a range, 
the spin magnetization varies within 2%: it is not likely 
to be responsible for the MAE increase. On the other 
hand the anisotropy of orbital moment strongly corre- 
lates with MAE and is probably responsible for its en- 
hancement as tetragonality increases. Orbital magnetic 
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FIG. 7: (Color online) K and the A as functions of c/a 
in Fei6N2. The ideal crystal structure without strain has 
c/a=l.l. For each c/a, the atomic positions are relaxed with 
volume being conserved. 



moments can be measured more precisely, so new XMCD 
type of experiments for this system are desirable. 



IV. CONCLUSION 

In this study of intrinsic magnetic properties of Fe 16 N 2 , 
our LDA results for magnetization agree with previously 
reported values while QSGTU increases magnetization by 
9%. This enhancement is largely due to on-site exchange 
splitting between the d minority and majority states - an 
effect seen in m any other magnetic systems such as NiO 
and MnA a ' °l In Fei6N2 in particular, we find no evi- 
dence of localized states or correlations not already found 
in Fe. Taken together all of those factors we expect that 
the QSGIU prediction for M is not far from what should 
be observed in the ideal Fei6N 2 compound. We find no 
evidence of charge transfer between different Fe sites as 
proposed elsewhere. Thus, the theoretical magnetization 
predicted for FeigN2 does not exceed the maximum on 
Slater-Pauling curve (~2.5/is) and is smaller than corre- 
sponding maximum of magnetization observed in Fe-Co 
alloys, which may still be considered as a record holder 
among d atomic magnets. 

LDA calculations predict Tq significantly larger than 
the experimental value; the QSGJU result is even larger. 
We assume that FeigN2 will have a higher Tq if one can 
find a way to stabilize it. Effects of doping by various 
elements on M and Tq were studied in the LMTO-CPA 
approximation. Various dopants affect M and Tq differ- 
ently; but unfortunately no dopants we considered en- 
hanced M or T c . 

A uniaxial magnetocrystalline anisotropy 

iiT=103xl0 5 erg/cni 3 was calculated in the LDA 
with the theoretically optimized crystal structure. K is 



strongly correlated with the atomic magnetic anisotropy 
energy due to spin-orbit coupling only. We found it 
can be increased by increasing c/a or by adding small 
amount of Co or Ti atoms on 4e or Ad sites. 

Fei6N2 is one of the more promising candidates for per- 
manent magnets that do not contain rare-earth elements. 
We believe that there is room for improvement and we 
studied several possible routes to obtain better proper- 
ties. A further investigation on increasing the thermal 
stability and/or changing crystal structure tetragonality 
is desired. 
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